Preprint typeset using I^TjrjX style cmulatcapj v. 5/2/11 



EVIDENCE FOR A WIDE RANGE OF UV OBSCURATION IN Z ~ 2 DUSTY GALAXIES FROM THE 

GOODS-HERSCHEL SURVEY* 

Kyle Penner 1 , Mark Dickinson 2 , Alexandra Pope 3 , Arjun Dey 2 , Benjamin Magnelli 4 , Maurilio Pannella 5 , 
Bruno Altieri 6 , Herve Aussel 5 , Veronique Buat 7 , Shane Bussmann 8 , Vassilis Charmandaris 9 ' 10 ' 11 , Daniela Coia 6 , 
Emanuele Daddi 5 , Helmut Dannerbauer 12 , David Elbaz 5 , Ho Seong Hwang 8 , Jeyhan Kartaltepe 2 , Lihwai Lin 1 ', 
Georgios Magdis 14 , Glenn Morrison 15 ' 10 , Paola Popesso 4 , Douglas Scott 17 , and Ivan Valtchanov 6 



Or 
D 

c/3 



ABSTRACT 

Dusty galaxies at z ~ 2 span a wide range of relative brightness between rest-frame mid-infrared 
(8 /im) and ultraviolet wavelengths. We attempt to determine the physical mechanism responsible 
for this diversity. Dust-obscured galaxies (DOGs), which have rest-frame mid-IR to UV flux density 
ratios > 1000, might be abnormally bright in the mid-IR, perhaps due to prominent AGN and/or 
PAH emission, or abnormally faint in the UV. We use far-infrared data from the GOODS- Herschel 
survey to show that most DOGs with 10 12 L Q < Lir < 10 13 L Q are not abnormally bright in the 
mid-IR when compared to other dusty galaxies with similar IR (8-1000 fim) luminosities. We observe 
a relation between the median IR to UV luminosity ratios and the median UV continuum power-law 
indices for these galaxies, and we find that only 24% have specific star formation rates which indicate 
the dominance of compact star-forming regions. This circumstantial evidence supports the idea that 
the UV- and IR-emitting regions in these galaxies are spatially coincident, which implies a connection 
between the abnormal UV faintness of DOGs and dust obscuration. We conclude that the range in 
rest-frame mid-IR to UV flux density ratios spanned by dusty galaxies at z ~ 2 is due to differing 
amounts of UV obscuration. Of galaxies with these IR luminosities, DOGs are the most obscured. We 
attribute differences in UV obscuration to either: 1) differences in the degree of alignment between 
the spatial distributions of dust and massive stars, or 2) differences in the total dust content. 
Subject headings: Galaxies: high-redshift — Infrared: galaxies — Galaxies: ISM 
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1. INTRODUCTION 

At z ~ 2, a larg e fraction of all high mass stars form 
in du sty galaxies (|Chapman et al.l [2005: Magnelli et al.1 
1201 If) . Most of the intrinsic UV emission from newly 
formed stars in these galaxies is obscured, or absorbed 
by dust grains that subsequently heat up and radiate at 
IR wavelengths. The IR luminosity resulting from this 
obscuration is usually much greater than the emergent 
UV luminosity. For galaxies in which the intrinsic UV 
emission from newly formed stars is less obscured, the 
IR luminosity is still greater th an the emergent U V lu- 
minosity, but to a lesser degree (jReddv et al ll2012t ) . The 
relation between the IR and emergent UV emission from 
a z ~ 2 galaxy depends on the interplay between star 
formation and dust obscuration. 

One of the many ways to select dusty galaxies at z ~ 2, 
without redshift determinations from spectroscopy, is to 
use the r atio of observed 24 to 0.65 /zm (i?- band) flux 
densities (TDev et al.l 12051 iFiore etaLl 120081 ). Sources 
satisfying S^/So^s ^ 1000 have been termed "dust- 
obscured galaxies", or DOGs; their redshift distribu- 
tion is approxi mately a Gaussia n that peaks at z = 2 
with cr 2 = 0.5 (jDev et al.l [20081) . In the redshift range 
1.5 < z < 2.5, 0.65 /im observations are sensitive to rest- 
frame UV emission from newly formed massive stars, and 
24 /im observations are sensitive to mid-IR emission from 
hot dust and polycyclic aromatic hydrocarbons (PAHs). 
The DOG criterion is thus unique in that it selects galax- 
ies in a specific redshift range, that also exhibit extreme 
ratios between their rest-frame mid-IR and UV flux den- 
sities. We have yet to understand the physical mecha- 
nism driving the span of ratios exhibited by dusty galax- 
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ies at z ~ 2. 

The IR luminosities of DOGs with Lir > 10 13 L Q are 
dominated by emission from active galac tic nuclei (AGN; 
iDev et aLlfeOOa iBussmann et alJl2009b[ ). The dominant 
sources of the IR l uminosities of less lumi nous DOGs is 
a topi c of debate. iFiore et all (|2008D and iTreister et al] 
( 2009) conclude that the IR luminosities of many DOGs 
with 10 12 Lfn < L T r < 10 13 Lq originate from AGN, 
while IPope et al l (pOOSf) conclude that many such DOGs 
are powered by newly formed stars. 

In this paper, we pose the question "What makes a 
DOG a DOG?" The primary goal of our study is deter- 
mining why DOGs have an extreme ratio between their 
rest-frame mid-IR and UV flux densities when compared 
to other dusty galaxies. Unfortunately, the simple and 
singular selection criterion cannot distinguish between a 
DOG that is: 

• abnormally bright at rest-frame 8 /im for its far- 
IR flux density, indicating its mid-IR luminosity 
may be dominated by AGN emission, or abnor- 
mally strong emission from polycyclic aromatic hy- 
drocarbons (PAHs); 

• or, abnormally faint in the rest-frame UV for its 
optical flux density, indicating that dust more com- 
pletely obscures the newly formed stars in the 
galaxy. 

We use Herschel (jPilbratt et al.ll2~010|) data in the Great 
Observ atories Origins Dee p Survey-North (GOODS-N) 
region (jElbaz et all 1201 ID to show that, on average, 
DOGs with 10 12 L < L m < 10 13 L© are not abnor- 
mally bright at 8 /im, but are more UV faint than other 
galaxies with similar IR luminosities. The ratio between 
rest-frame IR and UV flux densities is set by the amount 
of obscuration, which can vary with either: 1) the degree 
of alignment between the spatial distributions of dust 
and massive stars, or 2) the total dust content. 

This paper is organized as follows. We present the data 
and sample selection in ^ in §5\ we show the results. 
We discuss the implications of these results in gl and 
conclude in £j5j We assume a cosmology with Hq = 70 
km s" 1 Mpc - 1 , fl m = 0.3, and fl A = 0.7. 

2. DATA 

2.1. Measured quantities 

Our study uses observations of the GOODS-N region, 
which is roughly 10 arcmin x 16.5 arcmin in extent. We 
cull the sample of DOGs from a catalog of 24 /im sources 
produced for the Spitzer /MI PS survey of the GO ODS- 
N region (M. Dickinson, PI; iMagnelli etaLl 120111 ). A 
24 /im source is defined as a > 3tr flux density measure- 
ment from PSF fitting to Spitzer /IRAC 3.6 /jm source 
priors. The catalog is 99% complete at S24 > 50 /zJy, 
and contains 1603 sources. 

The 2.2 /im (i4T s -band) image we use to identify coun- 
terparts for the 24 /im sources comes from observations 
using the Canada- Fr ance-Hawaii Telesco pe (CFHT) . The 
data are pr esented in Wang et al.l (|2010t ): we use our own 
reductions (|Lin et alj|2012j ). The 0.65 /mi (i?-band) Sub- 
aru image we use t o define the DOG sample comes from 
iCapak et all (|2004D . The 5a depth of the 2.2 /im image 
is ~0.60 /Jy (24.5 AB mag); the 3er depth of the 0.65 /im 
image is ^0.05 /iJy (27.2 AB mag). 



To extract flux densiti es, we follow a mo dified version 
of the p rocedure used by IPope et al.l (|2008l) . Using SEx- 
tractor (|Bertin fc Arnoutsl 119961 ). we place 3" diameter 
apertures at the positions of sources detected (> 5a) in 
the 2.2 /im image. If the 2.2 /im flux density is detected 
with S/N > 5a but the 0.65 /im flux density is not de- 
tected with S/N < 3cr, we use a 3a limit for the latter 
flux density. 

To determine rest-frame UV continuum power-law in- 
dices, we extract flux densities at 0.45, 0.55, 0.80, and 
0.90 /i m (the B-, V-, J- , and z-bands) from Subaru im- 
ages (jCanak et alJl200l . using the same procedure. We 
use the 3.6, 4.5, 5.8, and 8 /im flux densities already as- 
sociated with the 24 /im sources to determine whether or 
not their spectral energy distributions (SEDs) at these 
wavelengths behave as power laws; these flux densities 
come from a catalog produced for the Spitzer /IRAC sur- 
vey of the GOODS-N region, and will be included in cat- 
alogs accompanying the GOODS- Herschel public data 
release. 

For the optical/near- IR photometry, we calculate aper- 
ture corrections, defined as the ratios of total flux density 
to flux density in a 3" diameter aperture for point sources 
(non-saturated stars). We take the SExtractor parame- 
ter FLUX_AUT0 as the total flux density. The corrections 
are factors of 1.086, 1.225, 1.247, 1.057, 1.086, and 1.057 
at 0.45, 0.55, 0.65, 0.80, 0.90, and 2.2 /im, respectively. 
To maintain the signal-to-noise ratios given by SExtrac- 
tor, both the flux densities and their uncertainties are 
multiplied by these factors. 

We associate each 24 /im source with a 2.2 /im source 
and its extracted optical flux densities if the 2.2 /im 
source is a unique match within 0.76" of the position 
of the 3.6 /im prior. The match radius is chosen by max- 
imizing the number of unique matches while minimiz- 
ing the number of multiple matches. Of the 1603 24 /im 
sources, 87 either do not have a > 5a 2.2 /im counter- 
part within the match radius (85 of 87) or have multiple 
counterparts (2 of 87). 

The far-IR flux densities come f rom a catalog prod uced 
for the GOODS- Herschel survey (|Elbaz et alJl2Qlll ). We 
only us e 100 and 160 / i m flu x densities measured with 
PACS (jPoglitsch et ~aLl I2010D . and 250 urn fl ux densi- 
ties measured with SPIRE (jGriffin et aT1l2010D . to avoid 



the complications of measuring flux densities for 24 /im 
sources in the 350 and 500 /im SPIRE images that are 
affected by severe source confusion. We consider a > 3a 
measurement at 100 or 160 /im to be a detection; at 
250/im, we require a > 5a measurement. 

We impose additional constraints on the 250 /im flux 
densities (and 5 a limits), similar to those defining the 
"clean index" (jHwang et al.l 120101 : lElbaz et al.l I2010L 
l20ll . The 250 /im flux densities of clean sources (and 
5a limits of clean non-detections) should not be affected 
by severe source confusion. We do not impose addi- 
tional constraints on the 100 and 160 /im flux densities 
because the 100 /im images are not deep enough to be 
affected by source confusion, and the 160 /im images are 
only deep enough to be moderately affected. The 3cr 
depths of the 100 and 160 /im images are ^1100 /iJy and 
~2700 /iJy, respectively; the 5a depth of the 250 /im im- 
age is ^5700 /iJy. 

We attempt to match each 24 /im source to a source 
with a spectroscopic redshift from the catalogs of 
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iBarger et al.l (|2008l which includes redshifts compiled 
from the literature) and Stern et al. (in preparation). We 
find spectroscopic redshifts for 910 (57%) of the 24 /im 
sources. If no coincident sources with spectroscopic red- 
shifts are found, we resort to the photometric redshift 
source catalog of Pannella et al. (in preparation) to find 
a source match. We exclude sources with photometric 
redshifts derived from ill-fitting templates (those with 
reduced x 2 > 2). For an additional 510 (32%) of the 
24 /im sources we have photometric redshift estimates. 
There are no redshift estimates for 183 (11%) of the 
24 /im sources, and these sources are excluded from our 
samples. 

2.2. Samples 

Using the multi-wavelength information and redshifts, 
we define 2 samples from the superset of all 24 /im sources 
with S24 > 50 /iJy: 

• DOG sample: All 24 /im sources with 2.2 /im coun- 
terparts and 5 , 2 4 /5'o.65 > 986 (Fig. [TJ. The 24 /urn 
flux density of the faintest DOG is 53 /iJy, justify- 
ing our 5*24 > 50/iJy cut for the control sample. 
The limiting quantit y is the depth of th e 0.65 /im 
image (3<r = 0.05 /iJy iCapak et aT1l2004D . The red- 
shift distribution of our sample of DOGs is shown 
in Fig. [2j also motivating our 1.5 < z < 2.5 cut 
for the control sample. Six (of 61; 10%) DOGs 
have spectroscopic redshifts. In the following anal- 
ysis, we include only DOGs with 1.5 < z < 2.5. 
Our conclusions do not change if we include DOGs 
without redshift estimates in the sample. 

• Control sample: All S24 > 50 /iJy 24 /im sources 
with 2.2 /im counterparts that are at 1.5 < z < 2.5, 
and that do not satisfy the DOG selection criterion. 
Seventy four (of 268; 28%) control galaxies have 
spectroscopic redshifts. 

For each sample, Table Q] characterizes the subset of 
sources with flux densities detected at 0.65, 100, 160, 
and 250 /im. More than 70% of these galaxies are un- 
detected in optical spectra (or are unobserved) because 
their observed- frame optical flux densities arc so faint. 

Our sample con tains DOG s with fainter 24 /im emis- 
sion than does the iDev et aTl (|2008[ ) sample. Their sam- 
ple is selected from the shallower Spitzer /MIPS survey 
of the Bootes region; their limit is S24 > 300 /iJy. How- 
ever, the GOODS-N region is much smaller than the 
Bootes regi o n, so w e have few DOGs with S 2 4 > 300 /iJy. 
iPope etaU (|2008f ) also study a sample o f DOGs in 
GOOD S-N. The main differences between the lPope et al.l 
(2008) sample and ours are that: 1) they limit their sam- 
ple to S24 > 100 /iJy; and 2) they estimate a redshift for 
each DOG using IRAC and MIPS photometry, whereas 
we match DOGs to sources in a near-IR/optical catalog, 
with redshift estimates based on UBVRIzJK, 3.6, and 
4.5 /im photometry. 

The fractions of 24 /im sources at 1.5 < z < 2.5 that 
meet the DOG criterion increase with increasing 24 /im 
flux density (Fig. QJ. Of the sources with 624 < 100 /iJy 
4% are DO Gs. Of the sources wi th S 2 4 > 100 /iJy, 25% 
are DOGs. iRiguccini et al.l ()201lD find similar fractions; 
they also find that of their sources with S24 > 300 /iJy, 
60% are DOGs. 
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Figure 1. 24 fim flux density vs. 0.65 /im flux density for galaxies 
in the two samples. The DOG sample is defined by S24/S0 65 > 
986. 



100: 



<~\- I 
1 i 1 

! i.-" 



10. 



DOG sample 
24u.m sources 




0.0 



0.5 



1.0 



1.5 
Redshift 



2.0 



2.5 



3.0 



Figure 2. Redshift distributions for the DOG sample and all 
24 fim sources (before we impose redshift limits). For the analysis, 
we limit all samples to the redshift range 1.5 < z < 2.5 (the area 
between the vertical lines). 

2.3. Derived quantities 

Several quantities are useful in analyzing the relation 
between IR and emergent UV emission from galaxies. In 
this section, we detail how we estimate the total IR and 
UV luminosities, UV continuum power-law indices, star 
formation rates, and stellar masses for the galaxies in our 
samples. 

2.3.1. IR luminosities 

We estimate a total IR luminosity (8-1000 /im; Lir) for 
each gal axy with detected emi ssion at 100 /im. We red- 
shift the lCharv fc E lbaz ( 2001) template spectral energy 
distributions (SEDs) to the distance of each galaxy, find 
the SED that most closely matches the observed 100 /im 
flux density, and multiply the IR luminosity of that SED 
by the ratio between actual and predicted flux densities 
to get L m . 

We prefer this approach over estimating the IR lumi- 
nosity directly, by summing several far-IR flux densities, 
because the latter procedure requires detected emission 
at 160 and 250 /im. The GOODS- Herschel image at 
160 /im is moderately affected by blending due to source 
confusion, while the image at 250 /im is so deep that 
blending is problematic. The drawback to our chosen 
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Table 1 









Number detected 


Sample 


Number with 
1.5 < 2 < 2.5 


Median z Median 524 0.65 Aim 
/jjy > 3a 


100 Aim 160/im 250 Aim 
> 3ct > 3(T > 5a, clean 


DOGs 
Control 


61 

268 


2.1 161 47 (77%) 
2.0 102 268 (100%) 


29 (48%) 24 (39%) 9 (15%) 
81 (30%) 52 (19%) 15 (6%) 


Note. 


— All sources 


have > 5<r 2.2 fim and > 3cr 24 Aim flux density measurements. We also 



require S24 > 50 n Jy. 

approach is that all statements we make regarding IR 
luminosities assume that the low-redshift template SEDs 
accurately represe nt the SEDs of galaxies at 2 ~ 2. 
lElbaz et al.l (|2010f ) show that this assumption is valid 
when template matching is done to 100 A^m flux densi- 
ties. 

2.3.2. UV continuum power-law indices and luminosities 

For galaxies with UV emission from newly formed mas- 
sive stars, the UV continuum can be approximated as a 
power law with an index /3: 



(1) 



We use > 3tr flux densities at 0.45, 0.55, 0.65, 0.80, and 
0.90 fixa (i.e., in the B-, V-, R-, I-, and z-bands) to fit for 
(3 and the constant factor C for each galaxy which has 
an estimate of Lir. If only two flux densities have S/N 
> 3ct, we calculate f3 analytically. If only one flux density 
has S/N > 3cr, we use an upper limit for the 0.45 jum flux 
density to calculate a lower limit to f3. 

We estimate a UV luminosity XL\ at rest-frame 
0.16 /im for each galaxy using its redshift estimate and 
the power-law fit to the rest-frame UV flux densities. 
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Figure 3. Distributions of IR luminosities (derived from best-fit 
templates) for galaxies with detected emission at 100 Aim. The two 
distributions are statistically indistinguishable. 
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2.3.3. Stellar masses and star formation rates 

To estimate a stellar mass for each galaxy, we fit stellar 
population synthesis models to its UBVRIzJK, 3.6, and 
4.5 Attn flux densities (jDrorv et al.l I2004t IPannella et al] 
l2009bD . Full details are in §4.2 of lMullanev et all (12012D. 
We as sume the stellar initial mass function in ISalpeterl 
(1955) fro m 0.1 to 100 Mp, as w ell as the dust attenua- 
tion law in iCalzetti et al.l (120001) . 

We use the equation in IKennicuttl (|1998|) to calculate 
a star formation rate (SFR), based on the IR luminosity, 
for each galaxy with detected 100 fj,m emission. We make 
no correction for the emergent UV luminosity, since it is 
negligible for all galaxies in the DOG a nd co ntrol sam- 
ples (see §3]). In using the IKennicuttl Jl998) equation, 
we assume that the observed 100 ^m emission is due to 
star formation and not AGN activity, a nd that the star 
forma tion episode lasts for < 10 8 years. IMullanev et al.l 
(2012) find that the former assumption is valid for most 
AGN with detected X-ray emission at z < 3. 

3. RESULTS 

Seventy six percent of the z ~ 2 galaxies with detected 
100 /xm emission have 10 12 L < L m < 10 13 L Q (Fig. 
[3]). The distributions of IR luminosities for the DOGs 
and the control galaxies are statistically indistinguishable 
{jp = 0.20 that the two samples are drawn from the same 
parent population, using a K-S test). 
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Figure 4. Stellar mass vs. redshift for galaxies with detected 
emission at 100 fim. The distribution of stellar masses for DOGs 
is not strongly different from that for control galaxies, though 
some control galaxies have lower stellar masses than do the DOGs. 
Twenty four percent of the control sample has M* < 5 X 10 10 Mq, 
as opposed to 4% of the DOG sample. 

Fig. [4] shows the stellar masses for DOGs and control 
galaxies with detected emission at 100 /im. The distri- 
butions are not strongly different (p — 0.02, using a K- 
S test), though some control galaxies have lower stellar 
masses than do the DOGs. 17 (of 72; 24%) control galax- 
ies have stellar masses M* < 5 x 10 10 Mq, while only 1 
(of 25; 4%) DOG has a stellar mass below this threshold. 

The infrared SEDs of most low redshift dusty galaxies 
peak between rest-frame ~ 60 and 100 A^ni. The rest- 
frame 8 /im luminosity is only a fraction of the total IR 
luminosity in these galaxies. To address whether or not 
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Figure 5. S100/S24 vs. redshift for galaxies in the two samples. 
Galaxies with limits at 100 /im are not shown in this figure, but 
we do include them in our statistical tests. DOGs are statistically 
indistinguishable from the control galaxies in S100 /S24 ■ 

DOGs are abnormally bright at rest-frame 8 /im, we re- 
quire a comparison of the rest- frame far-IR and 8 /im flux 
densities between DOGs and the control galaxies. 

Fig. [5]shows that DOGs are statistically indistinguish- 
able from the control galaxies when looking at the ob- 
served flux density ratio S100/S24 (p = 0.35, using a 
K-S test) . In this figure, we do not show the galax- 
ies with limits for Sioo- We perform the Gehan and 
logrank tests, which are conceptually similar to the two 
sample K-S test, but also allow for the inclusion of galax- 
ies with limits at 100 /im in the two samples. We find 
no statistically significant difference between the DOGs 
and the control galaxies (p > 0.70 that both samples 
are drawn from the same parent population, from both 
tests). Our conclusions are the same using 5i6o/5*24 (not 
shown; p > 0.10 from both tests). If DOGs were abnor- 
mally luminous at rest-frame 8 /im for their far-IR lumi- 
nosities, then we would expect them to have low values 
of observed 6*100 / S24 and/or S'i6o/<S'24 compared to those 
of the control galaxies; they do not. However, 100 and 
160 /im observations are sensitive to rest-frame 33 and 
53 /im emission from galaxies at z = 2. These rest- frame 
wavelengths are still short of the presumed wavelength of 
the peak of the infrared SED; the rest- frame luminosities 
are still only a fraction of the total IR luminosity in these 
galaxies. 

250 /im observations are sensitive to rest-frame 83 /im 
emission from galaxies at z = 2. This rest- frame wave- 
length is generally close to the wavelength of the peak 
of the IR SED. Fig. [5] shows that DOGs are statistically 
indistinguishable from the control galaxies when looking 
at the observed flux density ratio S^so/S^. We reach 
the same conclusion when including galaxies with limits 
at 250 /im in statistical tests (p > 0.76 that both samples 
are drawn from the same parent population, from both 
a Gehan and logrank test). 

The DOG criterion does not select galaxies that are 
abnormally bright at rest-frame 8 /im for their far-IR 
flux densities. What makes a DOG a DOG must be 
that the galaxy's rest-frame UV emission is suppressed, 
compared to the UV emission from a control galaxy. 
Therefore we expect a clear separation between DOGs 
and control galaxies when looking at the ratios of rest- 
frame optical to UV flux densities. Indeed, Fig. [7] shows 



Figure 6. S2TM/S2A vs. redshift for galaxies in the two samples. 
Galaxies with limits at 250 fim are not shown in this figure, but we 
do include them in our statistical tests. As in Fig. \E\ DOGs span 
the same range of flux density ratios as do the control galaxies; this 
leads us to conclude that the DOG criterion does not select dusty 
galaxies that are abnormally bright at rest-frame 8 fim. 
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Figure 7. S2.2/ 'Sb.65 vs. redshift for galaxies in the two samples. 
A DOG is distinguishable from a random galaxy with detected 
24 /*m emission because it is abnormally faint in the rest-frame 
UV. 



that 92% of DOGs have an observed S 2 . 2/S0.65 > 20, 
while 78% of control galaxies have S2.2fS0.Q5 < 20. The 
Gehan and logrank tests report a statistical difference 
between the DOGs and control galaxies (p < 0.0001 that 
both samples are drawn from the same parent popula- 
tion, from both tests). Galaxies with S2. 2/ So. 65 > 20 
are also known as "extremely red objects", or EROs 
(lElston et alil98llHu & Ridgwavll99llGraham fc Devi 
119961 : iMcCarthvl 120041 ). We'note that even though a 
galaxy with detected 24 /im emission may meet the ERO 
criterion, it may not meet the DOG criterion (and vice 
versa, rarely). 

The statistical difference between the DOGs and con- 
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trol galaxies for S2.2/S0.65, and the lack thereof for 
S'ioo/5'24, is robust against photometric redshift errors. 
We offset photometric redshifts for the superset of all 
galaxies with detected 24 /im emission, assuming the off- 
sets obey a Gaussian distribution with ctaz/(i+*) = 0.1. 
We then select new DOG and control samples. The dif- 
ferences in S2.2/S0.65 are always statistically significant, 
and the differences in S100/S24 are very rarely statisti- 
cally significant (from K-S tests). 

4. DISCUSSION 

We have shown that the 100, 160, and 250 to 24 fim flux 
density ratios for DOGs with moderate IR luminosities 
(10 12 L Q < Lir < 10 13 Lq) are statistically indistin- 
guishable from those flux density ratios for galaxies with 
detected 24 /im emission that lie at similar redshifts and 
have similar IR luminosities, but that do not meet the 
DOG selection criterion. Most DOGs have higher 2.2 to 
0.65 fixsx flux density ratios than do the control galaxies. 
Thus it seems clear that DOGs occupy the tail of a dis- 
tribution of UV obscuration in IR-luminous galaxies at 
z~2. 

We select a sample of DOGs using deep images of 
the GOODS-N region. These DOGs have lower IR lu- 
minosities than do DOGs selected using sha llower im- 
ages of wider regions. For example, the iDev et ahl 
(2008) sample in the wide Bootes region contains many 
DOGs with Lir > 10 13 L©. Our sample contains very 
few galaxies with Lir > 10 13 L Q , and our conclusions 
regarding obscuration may not apply to DOGs with 
such high luminosities. At low redshift, AGN emis- 
sion dominates th e IR luminosity fo r most galaxi e s with 
Lir > 10 13 L dTran et al.ll200l . IDev et all (12008?) 
find that many DOGs at z ~ 2 with Lir > 10 13 L© 
have featureless SEDs from observed-frame 3.6 to 8 /im, 
indicating t he presence of AGN-heated dust. Sev- 
eral studies dTvler et al.l 120091: iBussmann et al.l l2009bt 
iMelbourne et al.l 120121) find that, for these DOGs, the 
SEDs at rest-frame mid- and far-IR wavelengths are 
similar to those of low redshift galaxies with IR lumi- 
nosities dominated by AGN emission, such as Markar- 
ian 231. Only 2 (of 59; 3%) DOGs, and 8 (of 253; 
3%) control galaxies, have SEDs that increase (in vS v ) 
from observed-frame 3.6 to 8 /im, and indeed both are 
among the DOGs with the highest Lir, (at 4.8 x 10 12 and 
9.5xl0 12 L©). Two galaxies are not enough to allow us to 
rule out the possibility that DOGs with Lir, > 10 13 L© 
differ in their 100 or 250 to 24 /im flux density ratios from 
control galaxies with similar luminosities. 

Most galaxies in our sample have 10 12 L© < Lir < 
10 13 L©. Whether their IR luminosities are dominated 
by emission from A GN o r newly formed stars i s not clear. 
iFiore et al.l ([2001 and iTreister et all ((20091 ) find that 
the average DOG in this Lir range has an X-ray spec- 
trum with a power-law index they interpret as indica- 
tive of heavily obscured AGN emission. An obscured 
AGN woul d presumably cont ribute to the IR luminos- 
ity as well. IPope et ahl (|2008D find the same average X- 
ray index for a similar sample of DOGs, but also find 
PAH emission in the mid-infrared spectra of the DOGs; 
they conclude that heavily obscur ed AGN em i ssion can- 
not coexist with PAH emission. IPope et al] £2008) at- 
tribute the X-ray emission of the DOGs to X-ray bina- 



ries rather than AGN. In our sample, 5 (of 61; 8%) DOGs 
and 31 (of 26 8; 12%) control galaxi es have detected X- 
ray emission ([Alexand er et al.N 2003). For galaxies with 
Lir < 10 13 L©, we do not find any evidence that AGN 
emission at rest-frame mid-IR wavelengths is more or less 
common in DOGs than in the control galaxies: the frac- 
tions of DOGs and control galaxies with increasing SEDs 
between observed-frame 3.6 /im and 8 /im are similar, as 
are the fractions with detected X-ray emission. No mat- 
ter what powers their IR luminosities, DOGs in this Lir 
range are abnormally UV faint, and are probably the 
most obscured galaxies. 

We have shown that the brightness of the emergent UV 
emission in z ~ 2 dusty galaxies can vary considerably 
from galaxy to galaxy, even though the emergent UV 
emission is only a fraction of the IR emission. In a dusty 
galaxy at either low or high redshift, most of the intrinsic 
UV emission from newly formed stars and/or AGN is ab- 
sorbed by dust grains, which emit at IR wavelengths. We 
now ask: "Are the UV-emitting regions that we see from 
a DOG or a control galaxy spatially coincident with the 
IR-emitting regions?" If yes, the variation in emergent 
UV emission is due to UV obscuration. If no, the emer- 
gent UV emission may come from stars in a "normal" 
galactic disk, whereas the newly formed stars ultimately 
responsible for the IR emission are completely obscured, 
perhaps in a compact, nuclear star- forming region. In 
this scenario, the variation in the rest-frame optical to 
UV flux density ratios might be due as much to differ- 
ences in the stellar populatio ns from galaxy to galaxy a s 
to differences in obscuration ( Charmandaris et al. 2004). 

We cannot yet efficiently spatially resolve the IR- 
emitting regions for large samples of high redshift dusty 
galaxies. Two pieces of circumstantial evidence support 
the statement that the UV- and IR-emitting regions in 
these dusty galaxies are spatially coincident: 1) a mi- 
nority of DOGs and control galaxies deviate from the 
median specific SFR for galaxies at z ~ 2; and 2) DOGs 
and control galaxies define a continuous relation between 
their median IR to UV luminosity ratios and median UV 
continuum power-law indices. 

The specific star formation rate of a galaxy is the ratio 
of its SFR to its stellar mass. Most star-forming galax- 
ies at the same redshift fall on a tight relation between 
SFR and stellar mass; this relation is referred to as the 
"main sequence" (iNoeske et al.l 120071: lElbaz et al] [20071: 
Dadd i et al] [20071: iPannella et al] l2009al) . lElbaz et al l 
(1201 ID find that SFR and stellar mass are directly pro- 
portional at all redshifts, and present the evolution of 
the median specifi c SFR for galaxies at < z < 2.5. 
lElbaz et al] ([20111 ) also isolate the minority of galaxies 
with specific SFRs much higher than the median val- 
ues ("starbursts"). Here, we define starbursts as galax- 
ies with specific SFRs higher than 3 times the median 
value at their redshift. A minority of DOGs (5 of 25; 
20%) and c ontrol galax i es (18 of 72; 25%) are starbursts 
(Fig. HJ}. lElbaz et all ([20 111) find that individual low 
redshift starbursts have compact UV-emitting regions; 
the correspondence between deviation from the median 
specific SFR and the compact size of the UV-emitting 
region also holds for the average z ~ 2 starburst. If this 
correspondence holds for individual z ~ 2 galaxies, few 
DOGs and control galaxies would have concentrated UV- 
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Figure 8. Specific star formation rate (SFR over stellar mass) 
vs. redshift for galaxies with detected 100 fira emission. The solid 
line is the media n specific SFR for star-forming galaxies, as a func- 
tion of redshift {Elbaz ct al. 201 lj). The dotted line is three times 
the median value, which we use as a threshold for identifying star- 
bursts. 20% of DOGs and 25% of control galaxies are starbursts. 
DOGs deviate from the median specific SFR no more frequently 
than do the control galaxies. 



Figure 9. L]r over 0.16/im luminosity (IRX) vs. UV continuum 
power-law index for galaxies with detected 100 /an emission. 
We only show errors on the UV power-law indices. Filled hour- 
glasses and circles show the median j3 values, in bins of IRX; the 
vertical bars show the range of IRX over which the j3 values are 
binned . The "average LBG" black asterisk is from IReddv et al.l 
(2012), and shows a mea n and its dispersion. T he solid line is the 
local IRX-/3 relation from lOverzier et al. (2011). 



emitting regions. These regions in the typical DOG or 
control galaxy would be more widely distributed, and it 
would be less plausible that the UV-emitting regions oc- 
cupy parts of the galaxy not occupied by IR-emitting 
regions. That radio-emitting regions in high redshift 
galax ies are widely distribut ed indirectly supports this 
idea (Ruiopa karn et al1l2011|) . 

For starburst galaxies at low redshift. | Meurer et al.l 
(|1999l and more recently lOverzier et al.ll2011[ ) find a re- 
lation between the ratio of IR to UV luminosities (de- 
noted IRX) and th e power-law index o f the SED in the 
UV (denoted (3). (jMeurer et al.l 119991 refer to low red- 
shift galaxies with bright UV-emitting regions as "star- 
bursts" ; the term does not necessarily refer to galax- 
ies that deviate from the main sequence.) This rela- 
tion is generally interpreted to mean that the IR emis- 
sion originates as UV emission from newly formed stars, 
which is partially absorbed by dust. The dependence 
between the emergent UV emission, the dust emission 
in the IR, and (3, in the local IRX-/3 relation is con- 
sistent wotli_the_depejid^ce_of dust absorption on /? 
in the ICalzetti et al.l (|1994 120001) dust attenuation law. 
In other words, we think that galaxies that lie on this 
relation are thought to have spatially coinciden t UV- 
and IR-emitting regions. IReddv et al.l (I20T0I [20121) show 
that most Lyman break galaxies (LBG) at z ~ 2 lie on 
the local IRX-/3 relation. Howeve r, both z ~ 2 LBGs 
and the low redshift starbursts in IMeurer et al.l (|1999f ) 
are much less dusty than DOGs. Some sub-classes of 
dusty galaxies deviate from the local relation, perhaps 
because their IR emission is unrelated to their emer- 
gent UV emission dGoldader et al.ll200l ICh apman et al.l 
l2005t iPapovich et alj|2006t iBauer et"al1l201il) . 

Fig. EH shows t he local IRX-/3 relation from 
lOverzier et al.l ()2011D and the relevant quantities for the 
z ~ 2 dusty galaxies in our samples. Individual galaxies, 
in both the DOG and control samples, are found on ei- 
ther side of the relation. Because DOGs are so faint in 
the UV, it is difficult to precisely determine the power- 
law indices of their UV SEDs; the large uncertainties 



on /3 values preclude us from reaching firm conclusions 
about individual galaxies. The median j3 values for the 
DOG and control samples tend to lie to the left of the 
local IRX-/3 relation; DOGs do not appear to deviate 
more strongly from the relation than do control galaxies. 
These z ~ 2 dusty galaxies define their own IRX-/3 rela- 
tion. As noted earlier, specific sub-populations of dusty 
galaxies deviate from the local relation, but the trend 
that we find between the averages for a broader z ~ 2 
population has not been discussed. Because the median 
(3 value increases with increasing IRX, from the LBGs 
to the DOGs, the IR emission is not completely indepen- 
dent of the emergent UV emission in these galaxies. This 
reinforces our initial conclusion that z ~ 2 dusty galax- 
ies populate a continuum of UV obscuration, and that 
DOGs are simply the most heavily obscured galaxies. 

DOGs in the luminosity range spanned by our sample 
are more obscured than the control galaxies for a rea- 
son unrelated to their far-IR luminosity, because the Lir. 
distributions of the two samples are statistically indistin- 
guishable. In galaxies with dust heated by UV emission 
from newly formed stars, the amount of obscuration af- 
fecting that UV emission can vary due to either: 1) dif- 
ferences in the degree of alignment between the spatial 
distributions of dust and newly formed massive stars; or 
2 ) differences in the t otal dust content. 

IMeurer et "all ()1999D derive the local IRX-/3 relation by 
assuming a uniform screen of dust in the line of sight be- 
tween us and the n ewly formed stars (s ee figure 8 in 
ICalzetti et al.l 119941 for an illustration) . iGordon et al.1 
(|1997|) show that this assumption can be recovered when 
the dust is instead distributed in clumps around star- 
forming regions. Thus, DOGs might be more obscured 
than the control galaxies because DOGs have more dust 
clumps surrounding star-forming regions. 

Variations in UV obscuration might also be caused by 
differences in the physical properties responsible for the 
shape of the sub-millimeter SED. For example, increases 
in the mass of cold dust can result in greater UV obscu- 
ration. Increases in the cold dust mass can also result 
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in increases in sub-millimeter emission. Since the sub- 
millimeter luminosity is usually a small fraction of the 
total IR luminosity, the increase in total IR luminosity 
due to a higher cold dust mass would not be easily mea- 
surable between 100 and 250 (im. We cannot fully test 
this hypothesis without sub-millimeter luminosities for 
galaxies in our samples. The sensitivities of current sub- 
millimeter facilities limit comparisons of the dust content 
between DOGs and control gal axies to the most lumi- 
nous (Xtr > 10 13 L ©) galaxies (jBussmann et al.ll2009bt 
IMagdis et aLll20Tl . 

Currently available data do not allow us to discrim- 
inate between any hypotheses for the physical mecha- 
nisms responsible for either the differences in the de- 
gree of alignment between dust and stars or differences 
in the dust content. For instance, two plausible hy- 
potheses are that: 1) galaxy inclination might be re- 
sponsible for the patchiness of dust in the line-of-sight, 
and 2) many DOGs are merging galaxies, and some as- 
pect of the merging process creates differences in the 
degree of alignment. To falsify either we require high 
resolution rest-frame optical images; the Cosmic As- 
sembly Near-infrared Deep Extragalactic Legacy Survey 
(CANDELS) will ob tain such images us ing HST /WFC3 
in the near futu r e (iGrogin et al.l 120111). We note that 
iKartaltepe et al.l ()2012D and iSchawinski et~aTl (|2012t ) ex- 
amine rest-frame optical images of DOGs in GOODS-S 
and conclude that most are undistur bed disk galaxies. 
Furthermore, iNaravanan et al.l (|2010l ) simulate isolated 
galaxies that meet the DOG selection criterion, so DOGs 
are not necessarily merging galaxies on the basis of their 
extreme rest-frame mid-IR to UV flux density ratios. 

Finally, UV obscuration is affected by the prominence 
of the "dust bump" at rest-frame 2175 A in a galaxy's 
attenuation curve. At z ~ 2, rest-frame 2175 A is red- 
shifted into the i?-band filter. The presence of this fea- 
ture in the attenuation curves of high redshift galaxies 
has been controversial; the bump is found, though, in 
the att enuation curves for several samples of gal axies at 
z > 1 (INoll et al.ll200l IBuat et"aII[20TlL l20ll . How- 
ever, doubling the amplitude of the average bump leads 
to a reduction in the 0.65 /im flux density by a factor of 
1.3. This is not large enough to explain the spread of 
S2.2/S0S5 in Fig. El 

5. CONCLUSIONS 

1. We cull a sample of dust-obscured galaxies 
(DOGs), or galaxies with S'24/S'o 65 > 986 that are 
at 1.5 < z < 2.5, in the GOODS-N region. We 
use deep GOODS- Herschel data to compare the 
emission from DOGs with that from other z ~ 2 
galaxies with detected 24 emission. 

2. The DOGs in our sample span 10 12 L Q < Lir < 
10 13 L Q . DOGs and control galaxies, with detected 
100 /im emission, have similar distributions of IR 
luminosities and stellar masses. 

3. We compare the rest- frame far-IR and optical flux 
densities of DOGs with those of the control galax- 
ies. DOGs have extreme ratios of S24/S0.65 not 
because they are abnormally bright at rest-frame 
8 //m for their far-IR flux densities, but because 
they are abnormally faint in the rest-frame UV. 



4. DOGs and the control galaxies scatter around the 
median specific SFR established by z ~ 2 galax- 
ies falling on the "main sequence" ; 20% of DOGs 
have specific SFRs greater than 3 times the median 
value, and are thus starbursts. If UV-emitting re- 
gions in high-redshift starbursts are distributed as 
they are in low-redshift starbursts, few DOGs have 
compact UV-emitting regions. 

5. For both the DOG and control samples, the median 
rest-frame UV continuum power-law index (/3) at a 
given IR to UV luminosity ratio (IRX) is lower than 
the index predicted by the local IRX-/3 relation. 
DOGs do not appear to deviate more from this 
relation than do the control galaxies. Over more 
than a factor of 100 in IRX, the median j3 value 
for these galaxies systematically increases with in- 
creasing IRX. 

6. These pieces of evidence suggest that, for most of 
these galaxies, the UV- and IR-emitting regions are 
spatially coincident. Thus, the range in rest-frame 
mid-IR to UV flux density ratios spanned by dusty 
galaxies at z ~ 2 is due to differing amounts of UV 
obscuration. DOGs are the most heavily obscured 
galaxies. 

7. Differences in the amount of obscuration between 
DOGs and other dusty galaxies at z ~ 2 can be 
due to: 1) differences in the degree of alignment be- 
tween the spatial distributions of dust and massive 
stars, or 2) differences in the total dust content. 

Our samples do not have many galaxies with Lm > 
10 13 Li©. In DOGs with these IR luminosities, where 
AGN emission may be the dominant source of dust heat- 
ing, our conclusions about obscuration may not be valid. 

Further information about the nature of obscuration 
in these galaxies can come from measurements that spa- 
tially resolve UV- and IR-emitting regions. For instance, 
the Ha to H/3 line ratio varies with the amount of 
dust absorption, so its variati on across a galaxy traces 
the spatial dust distribution. iBrand et al.l (|2007f ) mea- 
sure the galaxy-wide Balmer decrement for several very 
luminous DOGs (those most likely t o have IR lumi- 
nositie s dominated by AGN emission): iMelbourne et all 
(|2011l ) present spatially resolved Ha maps, but H/3 is 
outside their spectral window. HST can resolve the 
rest-frame UV and optical emi s sion from these galaxie s 
(jBussmann et all l2009al 120111: IKartaltepe et all 120121) . 
and ALMA will resolve the IR emission and allow us 
to measure sub-millimeter luminosities we need to deter- 
mine cold dust mass. Studies using such data will clearly 
establish the degree to which the spatial distributions of 
dust and massive stars align. 

This work is based on observations made with Her- 
schel, a European Space Agency Cornerstone Mission 
with significant participation by NASA. Support for this 
work was provided by NASA through an award issued by 
JPL/Caltech. The research activities of M.D. and A.D. 
are supported by NO AO, which is operated by the As- 
sociation of Universities for Research in Astronomy un- 
der a cooperative agreement with the National Science 



UV obscuration in z ~ 2 dusty galaxies 



9 



Foundation. We are grateful to N. Drory for sharing 
the SED-fitting code we use to estimate galaxy stellar 
masses. 
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